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Abstract 

It is shown that the F4 rational and trigonometric integrable systems are exactly- 
solvable for arbitrary values of the coupling constants. Their spectra are found ex- 
plicitly while eigenfunctions by pure algebraic means. For both systems new variables 
are introduced in which the Hamiltonian has an algebraic form being also (block)- 
triangular. These variables are invariant with respect to the Weyl group of F4 root 
system and can be obtained by averaging over an orbit of the Weyl group. Alternative 
way of finding these variables exploiting a property of duality of the i<4 model is pre- 
sented. It is demonstrated that in these variables the Hamiltonian of each model can be 
expressed as a quadratic polynomial in the generators of some infinite-dimensional Lie 
algebra of differential operators in a finite-dimensional representation. Both Hamil- 
tonians preserve the same flag of spaces of polynomials and each subspace of the 
flag coincides with the finite-dimensional representation space of this algebra. Quasi- 
exactly-solvable generalization of the rational F4 model depending on two continuous 
and one discrete parameters is found. 
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1 Introduction 



The model was originally found in the Hamiltonian reduction method in ref.|jl| (see 
also the review Q). This model describes a quantum system in four-dimensional space 
with Hamiltonian which is associated with the root system R of the -F4 algebra |^]: 



1=1 j>i 
4 

I / \ I /zi + (-l)'^2 22 + (-l)'^3z3 + (-l)''4 24\ 

+ 5'2^^(2;i) +5 2^ ^(^ ^ ^ — ^-S^ — ^-J^ — ±j . 

j=l !^'S=0,1 

The arguments of a potential function v depend on the scalar products (a, z) = X]i=i c^j-^j) 
and summation goes over all positive roots a G = {ej, Cjibej, (ei±e2±e3±e4)/2, 1 < 
i,J <4}. 

In the rational model the function v takes a form 

<z) = 1 , (1.2) 



z2 



and the corresponding Hamiltonian becomes 
4 



oy(OP,r) 
«F4 - 2 
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" [zi + (-l)-2z2 + (-l)'^3 23 + (-1)^4^4]^ 



where a; is the harmonic oscillator frequency and 51,31/4 > —1/8 are the coupling con- 
stants. The configuration space is given by 

— 00 < Zi < Z2 < -^3 < -24 < 00 • (1-4) 

and it coincides with the Weyl chamber. 

In the case of the trigonometric i*4 model the oscillator term in ( |1.1| ) is absent, a; = 0, 



and 



v{z-l5) = -^, (1.5) 
sm pz 
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where /3 is a parameter. Hence, the Hamiltonian takes a form 



4-0 . -igai...^ 

§^1^^^+^^' ^sm^/? i-^+(-^)"-^+(-^)"-3+(-i)-^^4] ' 

with couphng constants 5,(71/4 > —1/8. The configuration space is given by the simplex 

vr 

—00 < Zi < Z2 < zs < Z4 < 00 , Zi — Z4 < — (1.7) 

y 

and it coincides with the Weyl alcove. When /? tends to zero the trigonometric Hamiltonian 
degenerates to the rational one at uj = 0. Both the rational and trigonometric F4 models 
are completely- integr able for arbitrary coupling constants g,gi. If g = the Hamiltonian 
(1.1) degenerates to the (rational or trigonometric) D4 model. 



Making a change of variables in (1.1) 

Zi^2 = Xi±X2 , Z3^4 = X3±X4 (1.8) 

we come to an equivalent Hamiltonian 
'^fM = ^Y.{-dl+4uj^xf)+2gY,Mxi-Xj;P)+v{xi + xf,p)} 

i=l j>i 
4 

+ 2gi^v{2xf,f3) + 2giY^ v [xi + i-irx2 + (-1^x3 + {-1^x4; l3] . (1.9) 

i=l u's=0,l 



Similarly to (1.1) this form can be associated with the dual root system of the F4 algebra 
i?Y = {2ei, Ci ± ej, (ei ± e2 ± 63 ± 64)} and we refer to the substitution ( pT^ ) as to 
the duality transformation. Due to relation v{2x; P) = (1/4)t;(x; 2/3) one can see that 
both in the rational and trigonometric cases this transformation interchanges the coupling 
constants gi ^ 2g. Besides that in the trigonometric case it rescales one part of the 
potential (this corresponds to changing in this part /3 — > 2/3 in potential functions v{x; (3)) 



preserving the other part. In particular, if = in (LI) the transformation ( |l.8|) converts 
( |1.1| ) into the Hamiltonian of D4 problem. 

It turns out that an analysis and formulas for the F4 trigonometric problem are much 
simpler in x-coordinates than in z-coordinates originally used. Similar phenomenon was 
observed in G2 model Q, where among two equivalent systems of relative coordinates (both 
of which were equally suitable for the rational model) the only one led to an algebraic form 
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of the trigonometric Hamiltonian ^ In what fohows we shall always use the form ( |1.9D of 
the Hamiltonian. 

In the present paper we demonstrate the exact solvability of the rational and trigono- 
metric -F4 models for general g,gi. The consideration uses a notion and a constructive 
criterion for exact solvability presented in This notion is based on the existence of 
a flag of functional spaces |^ with an explicit basis preserved by the Hamiltonian. A par- 
ticular criterion for exact solvability consists of checking whether the flag is related to 
finite-dimensional representation spaces of a Lie algebra of differential operators. If this 
criterion is fulfilled, then the Hamiltonian of the given system can be written in terms of 
the generators of this algebra which is called the hidden algebra of the system. In H it was 
shown that the eigenfunctions of the A^-body Calogero and Sutherland models g form 
an infinite flag of linear spaces of inhomogeneous polynomials in {N — 1) variables, which 
coincide to finite-dimensional representation spaces of the algebra gl{N), realized by first 
order differential operators in symmetric representation. The corresponding Hamiltonians 
were rewritten as quadratic polynomials in the generators of the maximal affine subalgebra 
of gl{N), and the coupling constants appear only in the coefficients of these polynomials. 
Recently, it was shown that this statement can be extended to all ABCD Olshanetsky- 
Perelomov rational and trigonometric integrable systems as well as to their SUSY gener- 
alizations which turned out to be associated to the hidden superalgebra gl{N\N — 1) (see 
Ref. 0). Later it was shown that for the 3-body Calogero-Sutherland models as well as 
BC2 models there exists a specific additional (second) hidden algebra which was called 
5*^^) C diff(2,R). This algebra turned out to be the hidden algebra of the G2 rational 
and trigonometric models ^ as well. Thus, the A2, BC2, G2 rational and trigonometric 
models are characterized by the same hidden algebra g^'^^ and their Hamiltonians can be 
written as non-linear combination of the g^"^^ generators. The flags which occurred in all 
above-mentioned models were always non-classical. Therefore their Hamiltonians were 
emerged in block-triangular form and the problem of spectra was reduced to a diagonal- 
ization of each separate block. However, a remarkable property of all above-mentioned 
models holds: a certain change of variables preserving the flag was sufficient to diagonalize 
all blocks simultaneously and, flnally, arrive at the pure triangular form. 

In the present paper we show that the Hamiltonians of both rational and trigonometric 
Fi models admit algebraic form, preserve the same inflnite non-classical flag of linear spaces 
of inhomogeneous polynomials and possess a hidden algebra which we can c diff(4,M). 
This algebra is an infinite-dimensional, finitely-generated algebra of differential operators 

^The form of the Hamiltonian is called algebraic, if exists, when the Hamiltonian can be represented as 
a Unear differential operator with polynomial coefficients 

sequence of linear spaces each one embedded into the next one, Vi (Z V2 ■ ■ ■ C Vn C . . . , forms an 
object called flag. A flag is called infinite flag (filtration) if a number of these spaces is infinite. A flag is 
called classical if AimVn = dimP„_i + 1. 
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possessing a finite-dimensional invariant subspace. It is worth to note that the D4 rational 
and trigonometric models, which play an important role in our analysis, are degenerations 
of i?C4 models from one side and -F4 models from another one. It implies that the D4 
rational and trigonometric models possess two different hidden algebras: gl{5) and f^'^^ as 
a consequence of these degenerations. 

The paper is organized as follows. In Section 2 the gl{5) Lie-algebraic form of the 
D4 rational and trigonometric models is studied. In Section 3 the rational F4 model is 
analyzed, its algebraic and Lie-algebraic forms are derived. In particular, a quasi-exactly- 
solvable generalization of the rational F4 model is found and studied. The trigonometric 
-F4 model is investigated in section 4. The variables providing the algebraic and the /^^^ 
Lie-algebraic forms of the Hamiltonian are found by averaging over an orbit of the Weyl 
group. We also discuss an alternative way to found such variables in connection with 
'dual' properties of the problem. Transition of the algebraic Hamiltonian from the block- 
triangular to pure triangular form based on introduction of new variables completes a 
demonstration of the exact solvability of the model. A realization of the algebra gl{5) 
in terms of first order differential operators acting on four-dimensional space is given in 
Appendix A. Appendix B is devoted to a description of the infinite-dimensional algebra of 
differential operators f^^\ admitting finite-dimensional representations in terms of inho- 
mogeneous polynomials in four variables. In Appendix C we give an explicit form of the 
variables leading to the algebraic form of the Hamiltonian in ^-representation. Finally, 
in Appendix D the explicit formulas for the first several eigenfunctions of the general F4 
model are presented. 



2 Algebraic and Lie-algebraic forms of the D4 rational and 
trigonometric models 

In this section we represent the Hamiltonians of the D4 rational and trigonometric models 
in an algebraic form, by making use of permutationally symmetric coordinates. The 
Hamiltonians can be written in terms of the generators of the gl{5) algebra ||^] and thus 
they have the gl{5) Lie-algebraic form. 



2.1 The rational D4 model 

The Hamiltonian of the D4 rational model is defined by (see 



1=1 



i<j 



+ 



(2.1) 
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where g = v{v — \) /2 > —1/8 is the couphng constant and oo is the harmonic osciUator 
frequency. The ground state eigenfunction of the Hamiltonian (|2.1| ) is given by 



exp(-*M) = (A_A+)''cxp|^-w^x?j , (2.2) 



with 



A± = \{{x^±x,) , (2.3) 

j<i 



and the ground state energy is 



£^0 = 4w(l + 6i^) . (2.4) 
As a first step towards an algebraic form of ( |2.1| ) let us make a gauge rotation 

/.g = -2(^W(x))-i(< - Eo)^t\x) . (2.5) 

Then, in order to code permutation symmetry and reflection symmetry Xi — > 

—Xi of the Hamiltonian [|, we take as new coordinates the elementary symmetric (Vieta) 
polynomials Si of the arguments x|, 

a I 2\ 2 I 2 I 2 I 2 

c/2\ 22, 22, 22, 22, 22, 22 

52 = S2[x)=Xi X2 +Xi Xs +Xi X4 + X2 X3 + X2 X4 + X3 X4 , 

a f 2\ 222, 222, 222, 222 

53 = ^3[X ) = Xi X2 X3 +Xi X2 X4 +Xi X3 X4 +X2 X3 X4 , 

54 = 5'4(x^) = Xi^ X2^ X3^ X4^ . (2-6) 

(r) 

In the s variables the Hamiltonian /i^^ becomes 

"S = E -^..^-l- + E(S, + A , (2.7) 

where 



■1 — ' — 7 -1 lyo^ 



k=l l>0 



These symmetries correspond to the group of automorphisms of the root space which in the D4 case 
is broader than the Weyl group. 
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A ■ ■ — A - ■ 

3 * — 'J ' 



^dxk'^J 2 ^ 
=2[1 + 2z.(4 - j)](5 - - So^js,- , (2.8) 

where so = 1 and Sj = at i > 4 or i < 0. 

Since the coefficients ( |2.8[ ) are polynomials in s, the operator ( |2.7] ) gives the algebraic 
form of the 1)4 rational Hamiltonian. It is worth to note that the coefficient matrix Aij 



makes sense of a flat space metric. An important feature of (2.8) is that the coefficients 
Aij and Bj are the second and the first order polynomials in s coordinates, respectively. It 
can be shown that it leads to two important conclusions: (i) the operator (2^) preserves 
the flag of spaces of polynomials 

piD,) ^ ^^Pi^P2^P3gP4| < pi +P2 +P3 +P4 < n) , (2.9) 
with the characteristic vector ^ 

/ = (1,1,1,1) , (2.10) 



and hence the operator ( |2.7| ) possesses infinitely many finite-dimensional invariant sub- 
spaces, and (ii) the operator (|2.7D has the gl{5) Lie-algebraic form Q, since it can be 
rewritten in terms of the gl{5) algebra generators but without raising generators (see 

(r) 

Appendix A). If = 0, u; = 0, the operator hj^^ becomes the fiat space Laplacian written 
in the gl{5) Lie-algebraic form. 

The operator ( p. 7] ) with coefficients ( |2.8D appears to be in pure triangular form with 
respect to the action on basis of monomials s^^ s^^ s^''' . Therefore, the spectrum of ( |2.7D , 

(r) 

hjJ^Lp = — 2e(/9, can be found explicitly and is equal to 

en = Mpi + 2P2 + 3p3 + 4p4) , (2.11) 

where n = 0, 1, . . . and pi are non-negative integers with a condition pi +P2+P?, +Pa = n. 
The spectrum does not depend on the coupling constant g, it is equidistant and corresponds 
to the spectrum of the harmonic oscillator. Degeneracy of the spectrum is related to the 
number of partitions of an integer number n to pi + p2 + p^ + Pa- The spectrum of the 
original rational D4 Hamiltonian (^]^) is = Eq + en- 

*A term 'characteristic vector' was proposed in [ p^ . It is a vector with components which are equal to 
the coefficients in front of p,; . 
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It is worth to mention that the boundaries of configuration space are determined by 
zeros of the ground state wave function (p.2|). In s- variables the boundary is an algebraic 
surface in four variables 

(A+A_)^ =Yl{x^ - x]f = 2bQsl - sl(192s3Si + 128s| - UAs2s\ + 27s\) 

i<j 



3„2\ 



+ 2S4(72S3S2 - Sspl - 4OS3S2SI + 9S3S2Si + 8S2 - 2S2S 

- 274 + 2sis^(9s2 - 2sl) - 52^1(452 -sl) = . 



(2.12) 



A simple relation between Jacobian and pre-exponential factor in the ground state wave 
function ( 2.12| ) exists 



det 



dsj 
dxk 



256 (A+A_) S4 . 



(2.13) 



Such a simplicity is of a general character because the Jacobian for the transformation from 
x's to the basic Weyl-invariant polynomials is equal (up to constant factor) to the product 
of linear functions vanishing on hyperplanes corresponding to roots (see, for instance, 



pil ), i.e. just to A+A_ from (2^). Since our s-variables (|2.6D differ from the D4 basic 
invariants in taking S4 = {xiX2X2X4^)'^ instead of X1X2X3X4, the extra factor S4 appears in 
the squared Jacobian. 



2.2 The trigonometric D4 model 

The Hamiltonian of the trigonometric model has the form Q 

1 1 



n 



(t) 
D4 



2^dx 

'1=1 



i<j 



+ 



sin {P{xi - Xj)) sin {P{xi + Xj)) 



(2.14) 



where g = z^(z^ — 1) /2 > —1/8 and /? is a parameter. When /? tends to zero the Hamiltonian 



(|2.14D coincides with (pl| at u; = 0. 

The ground state wave function is P, 0| 



with 



(A_(x,/3)A+(x,/3))^ = p-^^" Yl\sm^ Px, - sin^ Pa 

i<j 



A±(x, P) = p-^ Yl f^i^i ^ ^i) ' 

i<j 



(2.15) 



(2.16) 
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and the ground state energy equals 



(2.17) 



Using the same approach as in Section 2.1, we make a gauge rotation of ( 2.14 ) with 
the ground state eigenfunction as a gauge factor, h'j^l = -2^!-^{H^^^- Eq)^q. A straight- 
forward calculation leads to the operator 



15, - ^ 52 + 1/ ^ [cot {P{xi - Xj)) (a, - dj) + cot (/3(x, + Xj)) {di + dj)] (2.18) 

i=l i<j 



In order to code the permutation symmetry and reflection symmetry Xi - 

well as the periodicity of the Hamiltonian, we introduce new coordinates as the elementary 



symmetric polynomials of the trigonometric arguments (for definition, see eg. (2.6)) 

okix) = Skiyf) , (2.19) 



sin I3xi 



(2.20) 



which in the limit /3 ^ coincide with (2^). In these variables the Hamiltonian h^^^ 
becomes 



where 



«J=1 



4^(2/ + 1 + j - i)ai^i^i(jj+i 



l>0 
32 



4/5 [iaiaj + {i - j - 2)ai^iaj+i + {i - j - 4)fTj_20-j+2] for i < j 



A ■ ■ — A-- 



(2.21) 



B, + Cj{v) = 2[1 + 2^.(4 - j)](5 - j>j_i - A(5^j[aj + (7 - j>a,-+i] . 



(2.22) 



It is assumed that uo = 1, and fij = at z < or i > 4. The operator ( p.21[ ) with 
coefficients ( p.22| ) gives the algebraic form of the trigonometric Hamiltonian. It is 
worth to emphasize that Aij makes sense of a one-parametric fiat space metric, which is 



reduced to (2.8) at /? = 0. Similarly to what appeared in previously discussed D4 rational 
case the coefficients Aij, Bj are polynomials in ex's of second and first degree, respectively. 
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Hence, the Hamiltonian h^^^ can be written in terms of 5/(5) generators realized as first 
order differential operators 0| (see Appendix A), producing the gl{5) Lie-algebraic form 
of the -D4 trigonometric Hamiltonian. Furthermore, it can be easily verified that the 
gl{5) Lie-algebraic form of ( |2.14| ) does not contain the raising generators J^. Hence the 
operator ( p. 21 ) with coefficients ( |2.22| ) preserves the same flag of the spaces 

of polynomials but in a variables with the same characteristic vector ( |2.10| ). li v = 0, 
the operator /i^^^ becomes the flat space Laplacian written in gl{5) Lie-algebraic form 
depending on a single free parameter /? 0. 

The operator ( |2.21| ) with coefficients ( 2.22| ) appears to be in pure triangular form with 
respect to the action on basis of monomials ai^a2^a^^a'^^. Therefore the spectrum of 
( 2.21 ) , h^i^^ip = —2eip, can be found explicitly 

e„ =2/?2 [5piipi + 2p2 + 2p3 + 2p4) + 10p2(P2 + 2p3 + 2^4) + ISpsfe + 2^4) 

+20pl - 3{pi + 2p2 + 3p3 + 4p4) + 4z^(3pi + 5p2 + Gps + 6^4)] , (2.23) 

where n = 0, 1, . . . , andpi are non-negative integers with a condition {P1+P2+P3+P4) = n. 
The spectrum of the original trigonometric D4 Hamiltonian ( 2.14D is En = Eq + Cn- This 
completes a demonstration of the exact solvability of the trigonometric D/^ model. 

The boundaries of configuration space are determined by zeros of the ground state wave 
function ( 2.15| ). In the a variables the boundary is an algebraic surface in four-dimensional 
space 

(A_(x, /3)A+(x, /3)f = (3-^^ n (si^' " ^i"' f^^^if 



i<j 



16 cr4<T2 - 4 [afai + erg) - (80 a^aias - aja^ + 128 cr|) o-g 



+ 18(fTjfT3fT4 + 80-4(73 + 8aiaf + a'^(Ji)a2 - 21 0% 





21o\o\ - 



192ct| 0-1^3 + 256cr| 



(2.24) 



(cf. ( |2.12 )). We should emphasize that the algebraic surface ( p. 24 ) does not depend on the 
parameter /? and therefore equation (2.24) defines the same surface as in equation ( p. 12] ) 
but in the space of the fi-variables. It means that the configuration space of both rational 
and trigonometric D/^ problems is the same when is written in appropriate variables. 

There exists a trigonometric generalization of ( |2.13| ) giving a connection between Ja- 

®The operators ( ^.7| ) with coefficients (2.8) and (2.21) with coefficients (2.22) can be considered as two 



different deformations of Lapiacian which preserve the same flag (2.£). An interesting question is about 
existence of other non-trivial deformations of Laplacian preserving the same flag. 
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cobian and ( |2.24| ) 



det 



dxkj 



2 



256 (A+(/?)A_(/3))2 (1-/32^71 + /3V2-/?^<T3 + /3V4) . (2.25) 



In the limit /3 — > the equation (|2.25| ) becomes (|2.13|) . 

3 The rational F4 model 

We shah derive in this section the algebraic and Lie-algebraic forms of the rational F4 
models which lead to polynomial eigenfunctions. This fact together with the explicit 
calculation of the eigenvalues, exhibits the exact solvability of the model. 

3.1 Algebraic form 

The ground state of the rational F4 model 

«e = i E {-€. + 4.^4) + 2,j: (^-^ + ^-^) (3.1) 

i=l j>i ^ ■' J ' / 

+ 25iE:;r2+85i E 

can be written as 



4 

. i 

i=l ^ ^ + {.-\Y^X2 + {-ly^X^ + (-1)'^4X4] 



2 ' 



^«(x) = (A_A+)'^(AoA)'^exp^-a;E^.'l > (3-2) 

where g = v{v — l)/2, gi = ^(^ — 1), and 

4 

A± = \\{xi ±Xj) , 



i=l 



4 



Ao = 24n^^ ' 

A= n h + (-l)-'^^X2 + (-1)-'^='X3 + (-l)-"^X4] , (3.3) 
while the ground state energy is 

£^0 = 4a;(l + 6fi + 6iy) . (3.4) 
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The transformation ( |1.S| ) demonstrates the 'dual relation' between two parts of the wave 
function 

(A+A_)2(x) = (AoA)2(z) , 

(AoA)2(x) = (A+A_)2(z) . (3.5) 

General statement of the Hamiltonian reduction method is that any eigenfunction of 
( |3.1| ) can be written in a factorizable form as 

^{x) = ^1^\x)Pf,{x) , (3.6) 

where Pf^{x) is a polynomial in Xj's. The operator having these polynomials as eigen- 
functions can be obtained by gauge rotation of (|3.1|): 



4^-) = -2(^«(x))-i(<) - i?o)<^(x) . (3.7) 

In order to find variables supposedly giving an algebraic form to the original F4 Hamil- 
tonian (|3.lD we consider as a criterion of their choice the invariance with respect to the 



symmetries of the Hamiltonian 10|, i.e. under the group of automorphisms A of the -F4 
root space (in the F4 case this group coincides with the Weyl group W). The invariant 
polynomials of the lowest possible degrees generate the algebra of VK-invariant poly- 
nomials. These polynomials (denoted below as t'^^) can be found by averaging elementary 
polynomials (a,x)^" over some group orbit (we used the orbit Q generated by the root 
ei + 62, other orbits give algebraically related invariants): 

4""^ (x) = (a, xf' , a = 1, 3, 4, 6 . (3.8) 

The powers 2a = 2, 6, 8, 12 are the degrees of the group W. 

These polynomials written explicitly as polynomials in s (see (|2.(]| )) have a form 

if ^ = - 12s3 + 2s2Si + si , 

if ^ = 80s4 - 52s3Si + '^sl + sf , 

= - 346s?S3 + 20sl - 72OS4S2 + 1270s?S4 + 1652^ + 86s|s? 

- 122siS2S3 + 366s^ + sf . (3.9) 

The basis of ti^^ allows some non-linear transformations preserving the Weyl invariance 

(-1 ^ l-l ! 
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(n) 

^3 


->■ 


+ a3(tf) 


)^ 


An) 
E4 




+ a,{t?^ 




jn) 

'^6 









(3.10) 

where 03^4^6) ^4,6; cg, are any numbers. Supposedly, these transformations do not destroy 
an algebraic form of the operator hp^ if it appears in variables (3.8). 
Choosing 

*i = 4"^ , 

12 ^ 12^ ^ ^ 

80 ^ 30 ^ 3 48 ^ 1 ^ ' 

720*^ 288^^ ^ ^ 27^ ^ ^ ^ 1440^^ ^ 1080^^ ^ 

allows to simplify the variables to the form of polynomials of minimal possible degrees in 
s: 

ti = Sl, 

1 

^3 = S3 - ■^•si S2 


1 1 2 

t4 = Si - -S1S3 + — S2 

te = S4 S2 - ^ S2 - ^ s| + ^ Sl S2 S3 - ^ S? S4 ■ (3.11) 

The Hamiltonian (|3.7| ) written in the coordinates ( |3.11| ) has an algebraic form con- 
firming our expectation 



rj2 O 

hi'l=T.^-bg^^ +Y.iBa + Ca{f^,u)+Ca{u;))—, (3.12) 



where summation goes over a,b = 1, 3, 4, 6 and the coefficient functions are 
An=4ii ^13 = 12*3, 

Ai4 = 16t4, ^16 = 24*6, 

2 2 20 4 2 „ 

^33 ~ ~3 *1 ^3 + *1 ^4 J ^34 — — 2 *l ^4 + 8t6 ; 
^36 = 16 t4^ - 2 ) ^44 = — 4 t3 t4 — 2 ti tQ , 
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A^Q — —4titl — Gt^tQ , Aqq — —12 t^tl — 6 tit^tQ , 

Aba = Aab , (3.13) 



and 



Bi = 8 , = -tf , 



B^ = -At3, ^6 = -8*1*4- (3.14) 

The coefficients Aab have a meaning of elements of metric with upper indexes which 
corresponds to the flat space. Thus, the operator ( p. 12 ) with the coefficients Aab and Ba 



only {Ca{fJ',iy') = Ca{uj) = 0) is the flat space Laplace operator. 

Terms stemming from the potential part of the Hamiltonian are proportional to fi, 

Ci(/i, u) = A8{u + ^) , C3(/i, u) = -2 {2u + fi)ti, 

C4(/i,i/) = -12i/t3 , Cq{ii,v) = -12vtit4, , (3.15) 

(cf. ( p.l4D ) and also to to 

Ci{oj) = -4.ojti , C3(w) = -12w*3 , 

Ci{to) = -IQivU , C6(w) = -24cj*6 • (3.16) 

The operator ( |3.12|) with the coefficients Aab-, Ba, Ca(/u,z/), and Ca(w) represents the 
algebraic form of the i<4-rational model. It is easy to check that in the t-coordinates the 
Hamiltonian /ip^^ preserves a flag of polynomials V^^*\ given by 

^f^^ = (if if CCI < Pi + 2p3 + 2p4 + 3p6 < n) , (3.17) 
with the characteristic vector / 

/ = (1,2,2,3) , (3.18) 

(cf. ( p.lO[) ). This vector reminds the highest root among short ones in the root system 
of the algebra F4 written in the basis of simple root^. The flag ( |3.17| ) remains invariant 
under non- linear transformations ( 3.1C| ). The Hamiltonian also continues to be algebraic 



under these transformations. We should mention that the first study of the algebraic form 
of -F4 model was carried out in [^] using a set of variables with the only difference from 
our variables ( pl| ) in []. It was found the characteristic vector (1,2,3,5) in variance 

^°We thank Victor Kac for this remark. The similar statement holds for y4n-Calogero-Sutherland models 
where the flags are characterized by /* = (1, f , . . . , 1) and G2 models where f = (1, 2) 
^^It corresponds to ae — 0,ce, = 3/8, be, — 3/32, de, = 3/8. 
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to ( |3.18| ). However, since different variable is of the form ( |3.10| ), it preserves the same 
flag 'p(^-i) ( ^.17 ). It implies that the flag indicated in is not minimal in variance to 



the statement made in this article. 

We were able to find one-parametric algebra of differential operators for which there 
exists finite-dimensional irreducible representation marked by an integer value of the pa- 
rameter. Furthermore, the finite-dimensional representation spaces corresponding to dif- 
ferent integer values of this parameter form the infinite non- classical fiag (|3l7| ) see 



Appendix B). We call this algebra f^^\ Like the algebra g^"^^ (see [^), the algebra /^^^ is 
infinite-dimensional but finitely-generated. The F4 rational Hamiltonian in the algebraic 
form ( 3.12| ) can be rewritten in terms of the generators of this algebra. 



The operator ( 3.12| ) with coefficients ( ^.13D -( ^.16| ) appears to be in pure triangular 
form with respect to the action on monomials a^^a^a^^a^. One can find the spectrum 
of ( |3T^ ), h^^lif = -2e^, explicitly (cf. ( |2ll] )) 

e„ = 2 + 3p2 + + 6p4) , (3.19) 

where n = 0, 1, . . . , andpi are non- negative integers with a condition pi -|- 2p2 + 2p3 -|- 3p4 = 
re. The spectrum does not depend on the coupling constants g, gi, is equidistant and 
corresponds (with different degeneracy) to the spectrum of the harmonic oscillator as well 
as the rational D4 model. Finally, the energies of the original rational F4 Hamiltonian 
dp) are = Eq + en- 

Configuration space of the rational F4 model (p]^) is defined by zeros of the ground 



state eigenf unction, i.e. by zeros of the pre-exponential factor in (3^). These zeros also 
define boundaries of Weyl chamber (see |p). The squared pre-exponential factor can be 
written clS Si product of two factors. The first one 

{A+A^f = -l92tl + 256tl, (3.20) 

corresponds to the rational D4 model ( ^.ij) appearing at 51 = (cf. (|2.12D ). It looks much 
simpler than in [^]. The second factor 

( Ao A) ^ = - 3072 4 + 4096 4 - 2304 tj te - 432 4 + 3072 h h 
- 768 tl t4 te + 480 tj U - 192tf Iq - 8tf 4 
+ 16tj tl + 84 h U - \4 *6 , (3.21) 

corresponds to a case of the degenerate i<4 model, 5 = (as it was noted in Introduction 
it is equivalent to the -D4 model in dual variables). Thus, a boundary of the configuration 
space of the rational F4 model is confined by the algebraic surfaces ( 3.20| )-( |3.21| ) of the 



14 



third and seventh orders, correspondingly]^, being in total given by the algebraic surface 
of the tenth order. 



According to general theory (see ||TT[) the relation between Jacobian and the pre- 
exponential factor in the ground-state wave function ( \i.2\ ) is straightforward: 



det 



dxk 



^ 'A+A_)'(AoA)' . (3.22) 



4096 



In order to find eigenfunctions one can derive recurrence relations and then solve them 
out. For several first eigenfunctions it can be done explicitly (see examples in Appendix 
D). Similar to what was found previously for both the Calogero (^„-rational) model |^ ^ 
and the G2 rational model |13, ^ there exists among eigenfunctions a family $„(ti) which 



depends on a single variable ti. This fact was already used in ||T^ in order to construct 
quasi-exactly-solvable many-body generalizations of the Calogero model (for definition of 
quasi-exact-solvability, see for example |l5|). This family of eigenstates appears due to a 
fact that the coefficients ^n, Bi, Ci in ( |3.13| ) - ( p. 16 ) depend on the single variable ti 



only. Therefore, it is easy to verify that beside the flag 'p(^-i) (3. 17) the Hamiltonian h^^^ 



F4 



preserves another flag of polynomials V^^^ , deflned by the spaces 

Vi'Kh) = (tf |0<pi<n) . 

Since Vn^ C V^"^^ for any n, then the flag V^^^ C V^^'^^ . It leads to a degeneration of the 

(r) 

general spectral solution for the operator hp^ to an equation 

4ti$" - [4a; ti -8(6^ + 61/ + 1)]^>' = , (3.23) 
which can be solved explicitly, 

$W = 412^+12-+!) ^^ti) , En = -4un, (3.24) 

where L^"^ is a Laguerre polynomial in a standard notation, ra = 0, 1, 2, . . . . Existence of 
the flag V^^^ (which looks like a truly minimal flag) is a consequence of very degenerate 
nature of the F4 rational model. The Hamiltonian of the -F4 trigonometric model does 
not preserve this flag. Actually, there exist other 'degenerate' flags preserved by the F4 
rational Hamiltonian and thus other inflnite families of eigenstates depending on (^1,^3) 
or (ti,t3,t4) variables only, which can be found explicitly. They are also a consequence of 
degenerate nature of the i*4 rational model and a reader can easily investigate them. 

^■^In coordinates proposed in [^o[ these surfaces appear as of the eighth and tenth orders, correspondingly, 
leading to the total algebraic surface of the 18th order. We classify this choice of variables as non-minimal. 
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3.2 Quasi-exactly-solvable generalization of the F4 rational model 

(r) 

Above-mentioned remarkable property of h^^ allowing a family of eigenfunctions depend- 
ing on one variable leads to a possibility to construct a quasi-exactly-solvable (QES) 
generalization (for discussion see, for example, |jl5|) of the rational model. In order to 
do this we will use the same trick as was used in [0]. We look for a QES generalization 
of ( |1.1D of a form 

n'^f = ni:} + (to . (3.25) 

Let us make a gauge rotation ( 3.25| ) in the form ( p.?] ) and then require that the resulting 
operator possesses ti-depending family of eigenfunctions. It results in the equation 

j^iqes)^ = {Atidf^ - [Aujti - 8(6/i + 61/ + l)]ai - 2y('?^")}$ = -2e^ . (3.26) 



where the spectral parameter e is related to energy of the Hamiltonian (3.25) through 
En = Eo + En with Eq given by (U). 

Now one can pose a question under what condition on potential V^''^^\ the operator 
/jtees) jg Lie-algebraic. The problem is similar to which appeared in and the solution 
is the following: make a gauge rotation of h^'^^^^ in such a way that (i) to get rid off the 
potential V^'^^^\ and (ii) to obtain the sl{2) Lie-algebraic form 

= 4J° J" - 4u;J° + 2[n + 4(6^ + 6z^ + l)],r + 4a J+ - 47 J" . (3.27) 

Here 

J+ = rfSi - nti , jO = hdi - ^ , J' = di, (3.28) 

are the generators of the sl{2) algebra, and the potential y^*?*^*) should be 
yiqes) ^ ^2^3 _ 2aujtl + 2a[2n + 1 - 7 + 2(6/x + 61/ + l)]ti 
+ 27[7 + 1-2(6^ + 6. + !)] 
^1 

where the constant terms in potential are dropped off. Finally, we arrive in x-variables at 
the Hamiltonian (x^ = -^i^) 

4"-) = ly {-dl + Aco'x^) +2gy( --^ + , ^ 

F4 ^' ^\(Xi-XiY (Xi + Xi)^J 

i=l j>i ^ ■' J I ' 
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4 1 



i=l 



i/'s=0. 



^ [Xi + (-l)'^2a;2 + (-1)^3x3 + {-IY^Xa] 



2auj{-x.^Y + 2a[2n + 1 - 7 + 2(6^ + 6z^ + l)]x^ 



27[7 + 1 - 2(6/u + + 1)] 



(3.30) 



where we know (n + 1) eigenstates explicitly. Their eigenfunctions are of the form 



^W(x) = (A_A+r(AoA)'^(x2)7p„(x2) exp -c^x^ - ^(x^) 



2\2 



(3.31) 



where Pn is a polynomial of degree n. Hence we constructed the sl{2) QES deformation 



of the F4 rational model. If in (|3.3C| ) the parameter gi (and, hence, fi) vanishes, we arrive 
at the sl{2) QES generalization of the 1)4 rational model. The latter differs from the sl{5) 
QES deformation found in |16]. 

To conclude a discussion of the rational case, one can state that the rational F4 model 
admits the algebraic and also the f^^^ Lie- algebraic forms. Since the Lie-algebraic form 
of (|1|) with coefficients (gl3|)-(CT) contain no positive-grading generators, an infinite 
family of v, /^-depending polynomial eigenfunctions of ( p. 12 ) occurs. 



4 The trigonometric F4 model 
4.1 Algebraic form 

Let us consider the trigonometric F4 system now. Its Hamiltonian can be represented as 

= -]:i2^l+2gV^{x,(3) + ^V2{x,2P) , (4.1) 



where g = viv — l)/2, gi = ^(^ — 1), and 

1 



V^{x,P) = p^Y. 



4 



+ 



sin j3{xi — Xj) sin P{xi + xj) 



T^2(x,2/3) = 4/32^-^1— 
sm 2pxi 



+4/32 y: 



^ sin^ P [xi + (-1)^^2X2 + {-1Y-'X3 + i-lY^x^] 



(4.2) 



(4.3) 
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The ground state of the trigonometric F4 model has the form 

= iA+{x,(3)A.{x,(3)r{Aoix,2(3)A{x,2(3)r , (4.4) 



where 



A±{x,f3)=p-^Ylsmp{x^±Xj) , 

j<i 

Ao(x,2/3) =^-^JJsin2/3a;, , 

i 

A(a;,2/3) =/?-8 H sin (3 [xi + {-1Y^X2 + {-ir X3 + (-1^x4] . (4.5) 



u's=OA 



Here A±{x, (3), Ao{x, 2/3), A(x, 2/3) are the trigonometric anafogs of the factors appearing 
in the rational case (see ( |3.2| ) ) . In the limit /3 — > they coincide with those of the rational 
case, A±(x,0) = A-i-(x), Ao(x,0) = Ao(x), A(x,0) = A{x). The ground state energy of 
the Hamiltonian ( [4.1[) is given by (cf. ( |2.15 )) 



Eo = 4/3^(71/2 + 14^2 ^ ^g^^) ^ ^4_g) 



Guided by general theory 
Calogero-Sutherland models fs], 



and experience gained with previous studies of the 
G2 model Q and the rational F4 model (see Section 
3), let us check first whether there exists a family of factorized eigenfunctions of (4.1) of 
the type ^{x) = ^'q*^(x)Pp4(x), where the Pp^ are polynomials in some variables. If this 
is the case, there is a chance that, in accordance with the conjectures made in the 
trigonometric F4 model also possesses a hidden algebraic structure. In order to verify this 
we make first the gauge transformation of ( |4.1| ) with the ground state eigenfunction (4.4) 
as the gauge factor, 

4*] = -2{^^^\x)y\nil- E,){^^^\x)) ■ (4-7) 
Crucially important step is to find variables (if exist) which lead to an algebraic form 



of the Hamiltonian (4/7). In the rational case the relevant variables were the polynomials 
generating the algebra of the Weyl- invariant polynomials of Xj. This algebra allows 
the grading by means of homogeneous polynomials. In the trigonometric case we are going 
to deal with Weyl-invariant symmetric trigonometric polynomials in variables 

sm( I3xi) 

Vi = 



f3 ' 

which can be also represented by polynomials in 

o-fc(x) = Sk{yf) 
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where Sk are elementary symmetric polynomials. 

We impose an important requirement of a correspondence (the 'correspondence' princi- 
ple, see above) that in the limit (3^0 the new variables should coincide with the variables 
( 3.11 ) found for the rational -F4 model. Therefore, let us define the Weyl- invariant trigono- 
metric polynomials by averaging the elementary trigonometric polynomials over an orbit 
Q, generated by the root ei + 62, 



rf>(..« = Ef^^^^^^) . 0=1,3.4,6, (4.8) 



which reproduce the expression ( p.8| ) in the limit /? — > 0. These polynomials rjf^^ in terms 
of elementary symmetric polynomials ex's ( |2.19| ) - ( |2.20| ) look as 

{^) 2 2 

n = 0-1 - g/? 0-2 , 

r^^^ = - 12(73 + 2(T2(Ti +af + ^(36(74 - - 3(72(7?) 



= 80(74 - 52(73(71 + ^(7| + (7f + ^/^^ (24(74 (7i + 8(73(72 - 2(7^(71 - (72(7?) 

r^^) = - 720(74(72 1270(74(7? -h 366(7| - 122(73(72(7i - 346(73(7? + 20(7^ 86(7^(7? 16(72(7^ 
+ + (-864(74(73 - 28 56(74(72(71 + 432(74(7? 24(73(7| 1182(73(72(7? - 254(7^(7i 

9 

- 84(7^(7? - 9(72(7?) -/3^(864C7| + 952(74(7| - 864(74(72(7? - 538(73(7^(71 84C7| 

9 

32 

+ 72(7^(7? 15(7|(7f ) (216(74 (7|(7i + lAa^a'^ - lOajai - ba^af) 

(-288(74(71 + 8al + Walaj) - f^P^'^^n + ^P'^al (4.9) 

The basis of rjp^ allows some non-linear transformations preserving the Weyl invari- 
ance, which are more general than ( |3.10| ) 

rf^-rf)+(?„(r(^);/3) , a = l,3,4,6 , (4.10) 

where qa{T; (5) are polynomials in r's with /3-depending coefficients of dimensions (2 a) |^. 
^■^The dimension of Ta is equal to (2 a) while for /3 it is equal to ( — 1). 
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Following the same criterion as for rational case to have the variables in the form of 
polynomials of minimal possible degrees in CTq, we get 



n =0-1 — ^o"2 , 

7"3 = 0-3 - ^ 0-1 0-2 - 2/?^(cr4 - ;^cr|) , 

6b 

1 1 2 

r4 =£74 - - CJi 0-3 + — 0-2 , 

13 1 3 

rg = 0-4 0-2 - — O-^ - - cr| + - (Ti 0-2 0-3 - - (Ti C74 . (4.11) 

It is worth to show the relations between the variables rjp^ and t^. 



1 



256 

12r3 + rf - 8(3\llu - ^Tsn) + USP^^n + — /3^^6 



tP = 80t4 - 32r3ri + rf + 16/32(-28r4ri + 3T3Tf ) + ^p^ilOurf + 3r| - 28Te) 

2048 ^fi, , 4096 o 

+ ^/36(r6ri + r4r3) + ^/^'^l , 

Tj^^) = - 720T6 + 1000r4rf - dGTsrf + rf + 96t| + 8/?2(904t6Ti + 632r4r3 - 392r4rf 
- 96r|ri + IbT^rf) + QA[3^ {-2A8tqtI + 426r| - 360x4X3X1 + 35x4 r^^ + 15x|xf ) 

512 ^fi, 39 9 3x 4096 

+ — /3^(-144x6X3 + 56x6x3 - 552x|xi + 126x4X3x2 + 3x|) + ^/3^(-64x6X4 

+ 24X6X3X1 + 54x|x2 + 9X4X|) + ^/?10(5x6X4Xi + 3X|X3) 

131072 ^12/ 2 . 3n 
+ ^— /3''(t| + 6x|) , 

while the inverse algebraic relations ~ x^ in terms of xi^^ - do not exist. 

Quite surprisingly, the variables X4^6 in ( [1.11 ) contain no explicit dependence on /3. By 
construction in the limit /? — > the variables x^ coincide with the variables ta of ( 3.111 ). The 
polynomials ( 4.1l| ) are algebraically independent and generate a certain algebra S'^ {(3) 
of Weyl-invariant trigonometric polynomials. It is quite clear that {(3) is isomorphic 
to the algebra . It simply implies that the algebra has /3-parametric realization 

iHD. 

Finally, the gauge-rotated operator ( |4.7| ) in the coordinates ( [4.11 ) has an algebraic 



form 



o2 o 

4*] = + + ' a,^ = 1,3,4,6, (4.12) 



a, b a=l 
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where the coefficient functions are 



and 



An = 4ri - 4/3Vi2 - I/5V3 - , 

8 S2 
= 12r3 - -/3'(4rir3 + t^) - -^P^tit^ , 

^14 = 16r4 - y/5^rir4 - y/^'^'^e , 

^16 = 24 re - 20/3Vir6 - y/3V| , 
2 20 8 

^33 = -3 rg + y ri T4 - -/32 (18t| + t2 T4 + 12r6) , 

^34 = tI r4 + 8 T6 - ^,9^ (ri re + 12r3 r4) , 

^36 = 16 r| - 2 rf rg - ^/32(9r3 rg + n r|) , 
A44 = -4 r3 r4 - 2 n re - 24/3^1 , 
^46 = -4 ri r| - 6 rs re - 36/3^ r4re , 
ylee = -12r3r| - 6rir4re - 8/32(6r| + r|) , 

= , (4.13) 



Si = 8 - 8/3V1 , B, = -rl - I/3V3 - I/3V4 , 

S4 = -4r3 - y/3V4 , Se = -8rir4 - 56/3Ve . (4.14) 

The coefficients Aab have a meaning of elements of a metric with upper indexes which 
corresponds to the flat space for any value of the parameter (5. Hence, the operator ( |4.12 ) 
with the coefficients Aab and Ba (when Ca = 0) defines the flat space Laplacian in an 
algebraic form. At /? — > the expressions ( [4.13| )-( fl.l4 ) become ( p. 12 )^( ^3. 13) ). Terms 
stemming from the potential part of the Hamiltonian are proportional to /x, v 

Ci = 48(iv + /i) - 8/32 (5z^ + 6^)ri , C3 = -2{2v + fi)Tf - 16p'^{3iy + 5/u)r3 , 

C4 = -12i/r3 - 24/32(3zy + 4Ai)r4 , Cq = -12vtiTa - 48/3^(22^ + 3/x)r6 , (4.15) 
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(cf. i^). 

Eventually, the operator ( 4.12 ) with the coefficients ( [4.13D -( [l:.15 ) presents an algebraic 
form of the -F4 trigonometric model. It is straightforward to check that the operator (4.12) 
with the coefficients Aab, Ba, Ca from ( [4.13D -( p:.15D preserves the same flag of spaces of 
polynomials 'P^^*) (3^T^) as in the rational case (see ( [3.12| ))P^ Thus, it is evident that the 
Fi trigonometric Hamiltonian in the algebraic form ( |4.12| ) can be rewritten in terms of 
the generators of the /'^'^^-algebra. 



4.2 Duality 

Let us make now a short digression on a parallel between two equivalent representations 
of the Fi Hamiltonian related to two dual root systems mentioned in Introduction. These 
representations are connected by the 'duality' transformation D, eq.( |1.8| ), 

z = Vx , V'^ = -V , 

and we shall refer to two forms of Hamiltonian as to z- and ^-representations, correspond- 
ingly. The functional form of the Hamiltonian in these representations is remarkably 
similar (but do not coincide) because of similarity of two root systems. For the ratio- 
nal case the only change occurs in values of coupling constants g and gi while for the 
trigonometric case a rescaling of the parameter P has to be done for some terms in the 
potential. 

The dual transformation (|1.8| ) being applied to the Hamiltonian ( [4.1| ) provides the 
following correspondence: 



i=l i=l 

V,{x,P)^V2{z,f3) , 
V2{x,2f3)^4Vi{z,p) . (4.16) 

It implies that after the substitution (|l.8|) we arrive at the equivalent Hamiltonian in the 
Olshanetsky-Perelomov form [0] 



4r^(^) = -^E^' +M/^-im(^^/?) + ^^^V^^2(^,/3) , (4.17) 



i=l 



Note that according to the papers |l^ the rational and trigonometric F4 Hamiltonians preserve 
different minimal flags. Neither of them correspond to the flag P*^^*' (3.17). 
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with evident rescaling of the spectrum 

. (4.18) 

The factors of the wave function ( |4.4p corresponding to the potentials Vi and V2 also 
satisfy the 'duality' relations: 

(A+A_)(x,/3) = (AoA)(z,/3) , 

(AoA)(x,2/9) = (A+A_)(z,/3) . (4.19) 

At /i = (corresp. v = G) the F4 Hamiltonian reduces to the D4 Hamiltonians in x- 
(corresp. z-) representation. The relevant variables providing an algebraic form to the 
1)4 Hamiltonian are ak{x,l3) — symmetric polynomials in yf = sin^(/3xj)//3^ (see section 
Z.2, eqs. (pl9| ) - (g30D). The new variables Ta{x,(3) being polynomials in a{x,(3) give an 
algebraic form both to ly- and to ^-parts of the -F4 Hamiltonian. This means that one can 
expect existence of algebraic relations between dual variables Ta{z, (3) and Ta(x, f5). Indeed, 
there exist algebraic relations which express the variables Ta{z,P) through the Ta{x,l3): 

Ta{z,l3)=PaiT{x,(3)-(3) a = 1,3,4,6, (4.20) 

where Pa{T ', P) are polynomials in r's with /3-depending coefficients of the following form 

px(r;/3) = 2n-^/3Vf + ^/3S, 

P3(r; P) = -8X3 - + 4/3^ ( - 16r4 + ^rf + ^nrs) - ^-fp'r^r^ + , 

1 116 
P4(r; p) = 16t4 + + -rf - 4/?^ (2tiT4 + -rfrg) + -P^ {2tg + t|) , 

peir; P) = -Utq - 24t| - 8x1^4 - 2x^3 - 4 ^1 

3d 

- 4^2(8^3^^ _ 16T1T6 - 4rir| - rfr4 - ^r^^ T3) 

- 16/3^ (6r| + 2nr3r4 + ^rfre + ^rM) + |/36(r3r6 + ir|) . (4.21) 

The inverse relations for r(j;,/3)'s as functions of r(2;,/3)'s are not algebraic. However it is 
possible to express algebraically t(x, /3)'s through the r(z, /3/2)'s with the use of eqs.(^]20|): 

Taix, P) = Pa{T{V-^x, P); P) = Pa{T{z/2, /3); /3) . (4.22) 

The relations ( 4.22[ ) suggest another way of finding variables in which the F4 model 
acquires an algebraic form — one has to look for algebraic relations between ak = 
5fc(sin2(/?x,)//32) and = Sk{sm\Pzi/2)/P^): 

fa{<y;P)=ga{^\P) . (4.23) 
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Polynomial functions fa and ga can be used as new variables which give an algebraic 
form to the Hamiltonian in x and z representations, correspondingly. Historically we used 
namely relations (4.23) for finding relevant variables (4.11): 



Ta = /a(0-;/3) . 

Explicit formulas for ga{cr;f3) are given in Appendix C. One can see that the functions 
fa have much simpler form than the functions ga what gives a preference to the x- 
representation. 

4.3 Wave functions and energies of the trigonometric F4 model 



Configuration space of the trigonometric -F4 model is defined by zeros of the ground state 
eigenfunction ( |4.4| ), i.e. of the expression (A^-A- Aq A)^. This expression can be written 
as a product of two factors, (i): 



fA, A_ 



256rf - 192r«^ 



(4.24) 



which corresponds to the trigonometric D4 model appearing at gi = 0, and (ii): (Aq A)^, 
which corresponds to the degenerate {g = 0) i<4 model (cf. ( 3.2[)|) , (|3.21| )). In order to 
find the second factor we use the second relation (4.19). It tells that this factor is equal 
to A^ A^ written in z-variables and, thus, it can be expressed through r4^6(-2j/3). 



(AAo)2 = 256t4(z,/3)3 - 192t6(z,/3)2 . 



The use of eqs.(^), (|^ gives 



(4.25) 



A2 A2 



- ^^''-^6 + 8r°r3r4 + IQrfrl - 8tiV| - 192 rfrarg + 480T^r|r4 
768Tf r4T6 + 3072Tir3T| - 432r| - 2304T|r6 + 4096r| - 3072t|) 



3^3 



-2 2„ 



+ ^/3'(^/^6 - 3riV3T4 - erfrl + 3riV| + 96x1^3X5 - 246rfr|r4 
+ 288Ti^r4T6 - lUOrf tstI + 144Tir| - 1536ri r| + 960tit|t6 + 1536Tir| 

- 288T|r4 - 768T3r4r6) + — /3^( 246 rfr3r| + JGStiTsutq + 312Tir|r4 

- 2,24:tIt^tq - 2At^tI - 68x1^x4X6 + 30Ti^r|T4 - 192 r| re - Hrfrgre 

32 

- 672T2r| + 360tiV| - 2AtItI) - —(j^TlT^TQ - 3riV3r| - A^^tIt^t^ 

- Grf t| - 48Tf t| + 96r2T|r4 + 8x1^x3X4X6 + 312xi x|x| - 96xix|x6 - 48x| 

- 272x|x6 + 432x3x| - 384x3x|) + — /3^(xiV| + 24x^x3X4X6 - 72X1V3VI 

y 



.4^3 



.3 2„ 



24 



+ 180tiV|t6 - 648ri rgri - IUtitItq - 384Tir3r| + 288r|T4 + WbGririTe 



972t| + 768r4r|) - l^^/jio (rfrsrl + 6rir|r4r6 + 24rir4r| - 18r|r| 



54r3r|r6) + (3r|+8r6)r|, 
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(4.26) 



(cf. ( |3.21 )), which corresponds to the degenerate F4 model at 5 = 0. A boundary of the 
configuration space of the trigonometric F4 model is confined by the algebraic surfaces 
( 4.24[) , ( [4.26 ) of the third and eighth orders, correspondingly. It is quite surprising that in 
the rational case the corresponding surface defined by ( 3.21| ) is of the seventh order while 
in (|4.26| ) there exist the only two terms of the eighth order: (8/3)/3^ (r^rg — 3T-['T3r4). 

It is remarkable that in the F4 trigonometric case the relation between the Jacobian 
and the ground-state wave function has the same simple form as in the rational case: 



det 



dxk 



4096 



(A+A_)'(AoA)' 



(4.27) 



(cf-O)- 

Let us now proceed to finding the spectrum of the Hamiltonian ( [4.12| ). It is easy to 
see that the operator ( 4.12| ) with the coefficients ( [4.13| )-( [1.15| ) has a block triangular form 
in the r variables unlike pure triangular form which is needed in order to find eigenvalues. 
In general, in order to reduce this operator to pure triangular form it is necessary to 
diagonalize each block separately, doing it one by one. Surprisingly, in our particular 
problem it can be performed in full generality just by introducing unique set of new 
variables (!) 



pi = n 



1 



P3 = T3- ^(i "^rf 



Pi = T4 



Pe = tq 



(4.28) 



having the same dimension as in( [4.11 ). It is worth to note that this substitution becomes 
singular at /3 = 0, reflecting the non-existence of bound states for the rational 1)4 and F4 
models in absence of the harmonic oscillator term in potential. This coordinate transfor- 
mation is of the type ( [4.10D and hence leaves the flag ( 3.17 ) invariant. Thus, we arrive at 
a conclusion that among Weyl-invariant coordinate systems (of minimal dimension) there 
exists unique one which leads to pure triangular form of the Hamiltonian with respect to a 
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basis of monomials. We were unable to see a relation with variables introduced in |17, 18] 
in framework of a general study of trigonometric Hamiltonians based on root system which 
should guarantee triangular form. 

In the coordinates ( [4.28| ) the Hamiltonian ( [4.7] ) takes the form 



32 



128 



9 



3/03 - 2/3 ( -P4 + P3P1 



32 



+ ( 3/3- + 32p4 + 8p3/Oi - 24/3" - y/^Vej^^ip^ + 2 



3/3 ^pip4 



+ 24p6 - 4/?2 /9ip6 - 2p3P4 



9: 



pipe 



6/3-W3 + 16/32(p2 + -pg 



+ - 9/3-V1P3 + 16p6 - 32/3V3P4 + 6/3-2pip4 9; 



P3P4 



+ 
+ 



6/3-2 3^^^^ _ 4^^^g ^ 3g^2 _ 48/32 ^3pg 



5' 



P3P6 



- - ^/9"Vi (3/01P3 - I6P4) + 2pip6 - 4p3P4 - 24/3V1 

27 / \ 

+ — /3-V?P4 + I8/3-V1P6 + 3/3-2p4(3pi/>3 + 8/54] + 12p3/06 



5^: 



+ 64/3^4/06 



52 

P4P6 



^/3-ViP4 - 9/3-^6 ( P1P3 - 4p4 ) + 18p3P4 + 48/3>^ 



5^: 



+ I6/32 l + 3z^ + 5/i p3 



+ 24/3" l + 3z^ + 4/i p4 



-/3Vi(l- 




^P3 - 


)/34 


dpi - 



l + Qu + Qp] - 0''pi[l-Qp-bu] dp, - 3/3-2(l + 4i/ + 4^)pi 



^/3-^(l + 4z. + 4/i)pi + 12z.p3 

27 / 

— /3-V? - 9p-^i^pip3 + 6/3-2(5 + 61/ 



+ 6/i P4 + 



1 + 2z^ + 3/i /96 



5, 



P6 ' 



(4.29) 



and it seems simpler than the operator ( [4.12] ) with the coefficients ( [4.13D -( [4.15| ). It is 
easy to check that the operator ( [4.29 ) is indeed a triangular operator. It is evident that 
this operator can be rewritten in terms of the /^^^ -generators like it was for the operator 

iH). 

Using the representation ( [4.29| ) the energy levels of the Hamiltonian, /ip]<y5 = —2€ip, 
can be found explicitly and are given by 



en =4/32[pi(pi + 2p3 + 3p4 + 4p6) + 2p3(P3 + 2p4 + 3p6) + P4(3P4 + 8p6) 

+ 6pl + z^(5pi + 6p3 + 9p4 + 12p6) + 2^(3pi + bps + 6^4 + 9p6)] 



(4.30) 



where n = 0, 1, . . . , and quantum numbers pa are non- negative integers with a condition 
+ 2p3 + 2p4 + 3p6 = n. The spectrum of the original trigonometric F4 Hamiltonian (4.1) 
isEn = Eo + en (cf. (pD). 
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The explicit expressions for the first several eigenfunctions of ( 4.29| ) in p- variables 
are presented in Appendix D. It is worth to mention that the equations describing the 
boundary of the configuration space remain algebraic in p-variables and are given by 

(A+ A_) ^ = - 36p-^ pfpe + mp-'^plpl - 288(3-^ pip4Pe + 256^1 - 192pl , (4.31) 

and 

(AAo) ^ = - 72p-^ pi^ (3 p3^ + 8pe)+ 36/3"^ pi^ (24 pi ^6 + 9 pi^ps^ + 16 P4' 

- 144/?-2 pi (12 pi^pi + 8 P3P1P6 + 32 P6P4 + 3 psVi + 6 pi ^4^ 

3072 pe^ + 4096 p4^ + 864 pi^ps^ + 2880 p4Pi Vs^ + 2304 pepi ^3 
+ 768p3P4Vi + 7680p4PiV6 - 2304 P3V6 - 432 p3^ 

- 192/?2 f9pi2p3/04^ - 32piP62 - 3piP3^ - 2OP1P3V6 + 6p4P3^ 



+ 36 P4>1 + 16 P3P4P6J + 64/?^ (^60 P1P4P3' + 12 Pl>3'/06 - 2 P4>3' 

+ 9piV3^ - I6P4V6 - 32pi2p6^ + 160pip3P4P6) + ( - 54P4V3 

+ 108 P4 Ve/oi + 27 P3 V4 Vi + 34 P6P3^ + 48 pe ^3 + 6 P3^ 
256 / 

- — ( - 88 P3 V4/O6 - 64 P4P62 + 81 p4^ + 32 pipe ^3 - 24 p3 V4 

+ 12 P1P6P3') + 2048/?^V3/04' (p3' + 3 pe) + ^/3'' Pe' (3 P3' + 8 pe) . 

(4.32) 

Now the boundary of the configuration space of the trigonometric -F4 model is confined by 
the algebraic surfaces ( 4.31 )"( 4.32| ) of the fourth and sixth orders, correspondingly]^ unlike 
the r- variables where they were of the third and eighth orders, respectively (cf.( |4.31j )- 

dm)). 



5 Conclusion 

We have found that the general rational and trigonometric i*4 integrable models with 
two arbitrary coupling constants are exactly-solvable. After gauging away the ground state 
eigenfunction these models look very much alike when written in a certain Weyl-invariant 
variables. Their Hamiltonians preserve the same flag of the spaces of polynomials and both 

^^Being in total the algebraic surface of the tenth order. 
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models are characterized by the same hidden algebra - each Hamiltonian can be written as 
a non-linear combination of generators. It is an infinite-dimensional but finite-generated 
Lie algebra of the differential operators, which we call f^^^ algebra. It is quite interesting 
that D4 rational and trigonometric models possess two hidden algebras: gl{5) and /^^^. 
Similar situation takes place for A2, BC2 rational and trigonometric models as well as for 
G2 rational model. Their hidden algebras were gl{3) and g^'^\ respectively, however, G2 
trigonometric model had the only hidden algebra g^'^^ (see 

Present work complements previous studies where the algebraic and the Lie-algebraic 
forms as well as the corresponding flags were found for the rational and trigonometric 
Olshanetsky-Perelomov Hamiltonians of ABCD series (and their supersymmetric gen- 
eralizations) 01 and G2 model p. In order to conclude a study of the whole set 
of Olshanetsky-Perelomov integrable systems appearing in the Hamiltonian reduction 
method it is necessary to perform the same analysis for remaining £^6,7,8 integrable ra- 
tional and trigonometric models. We consider it as a challenging task for future. 



Our concrete consideration does not confirm some results presented in 17 1 for the 
rational and trigonometric models. In particular, we found that the minimal flags 
of spaces of polynomials preserved by the rational and trigonometric F4 Hamiltonians 
coincide, however, for both cases they differ from those given in [|l^, For trigonometric 
case our Weyl-invariant variables leading to pure triangular form of the £4 Hamiltonian are 
singular in (3. It reflects the fact that the triangular form for rational and trigonometric 
models has a different origin: for the rational case the diagonal matrix elements are 
proportional to ui while for trigonometric one to /3. The variables ( [4.28| ) look different 
from those presented in a general scheme |17, 18| formulated for all root systems. 



Acknowledgement 

J.C.L. and A.V.T. express their gratitude to M. Rosenbaum for the extensive discussions on early 
stage of the work (1997-1998). A.V.T. thanks N. Nekrasov for useful conversations and interest to 
the work. K.G.B. thanks ICN (UNAM) for kind hospitality where the collaboration was initiated, 
K.G.B. and J.C.L. also thanks LPT (Universite Paris-Sud) where the collaboration was continued, 
A.V.T. thanks ITAMP (Harvard University) for inspiring interaction with its members where the 
present work was almost finished and NSF grant for ITAMP at Harvard University for a partial 
support. The work was also supported in part by DGAPA grant No. IN120199, CONACyT grant 
25427-E, INTAS grant 00-00366, NATO grant PSTCLG 977275, and by Russian Fund of Basic 
Research 00-15-96786 and 01-02-17383. 



28 



A Representation of the algebra gl{5) 

The algebra gl{5) has a reahzation in terms of first order differential operators in four- dimensional 
(si, S2, S3, S4)-space: 

J- = d^, i = 1,2,3,4, 
J°k = Sidk , i. A: = 1,2,3,4 , 
J° = n - ^ Sidi , 

J+ = SiJ\ i = 1,2,3,4, (A.l) 

where n € C is a free parameter. In this realization a grading {ai\i = 1, 2, 3, 4) can be assigned to 
the generators (A.l) through their action on monomial 

7„PieP2 P3.P4 Pi+ai n2+a2 P3+a3„P4+a4 
u o-^ 02 ^3 La. o-]^ ^2 ^3 

Then the grading of J is defined as a four-component vector a = (ai, 02, as, 04). 

If n is a non-negative integer, the representation (A.l) becomes finite-dimensional and the 
corresponding representation space is a linear space of polynomials 

Vn = (sf sf I < (pi + P2 + P3 + P4) < n) . (A.2) 

For fixed n the algebra (A.l) acts on the space (A.2) irreducibly. As a function of n the spaces Vn 
possess a property that Vn C Vn+i for each n e Z+ and form an infinite fiag (filtration) 

U Vn=V. 



B The algebra 



We define the algebra /'^^ as an algebra of differential operators on which acts irreducibly on 
the space of inhomogeneous polynomials in four variables 

Vn = (sf 4" sf\0<pi+ 2p2 + 2p3 + 3p4 < n) , (B. 1) 

where neAf. Thus, /(^) C diff(C'*). 

The structure of the algebra /^^^ is the following. It contains three abelian subalgebras 
ij(fc), A; = 2, 3, 4 of first order differential operators 



>(2) 
■i 


= S\d2 , 


i = 0,l,2 , 


■i 


= S\d3 , 


i = 0,l,2 , 


■i 


= s{di , 


i = 0,1,2,3 



Rl"' =s\di , i = 0,l,2,3 , (B.2) 
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and 1 1-dimensional subalgebra B of first order differential operators 



4^) 


= di , 




= sidi , 


5^^) 


= 52^2 , 




= 53^2 






B^ 


= 53^3 , 




= 52^4 , 




= S3<94 


Bt' 


= S4(94 , 




= 818284 , 


r(4) 


= Si 53^4 • 







These algebras form a subalgebra B k (i?'^' 
second order differential operators 



®i?(4)) c /W. 



(B.3) 

Also there exists a set of the 



^(11) 


= S29?i , 


^(11) 
-^3 




7^(12) 
^4 


= S49i2 > 


7^(13) 
^4 


= 54913 ' 




= 84822 , 


^(22) 


= S1S49I2 , 


^(23) 
^4 


= S49I3 , 


7^(23) 
^14 


= S1S49I3 


^(33) 


= S49I3 , 


7^(33) 
^14 


= S1S49I3 , 


7^(14) 


— *2<^14 ^ 


7^(14) 
^23 


= S2S39i4 


7^(14) 
^33 


- s2f92 

— *3<^14 7 


7.(44) 
-^222 


— i2'^44 ; 


7^(44) 
^223 


— 52*3^^44 ^ 


7^(44) 
^233 


= S2S39I4 



(44) 
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■'3"44 J 



and third order differential operators 



4 — *4Wlll 



7.(222) _ 2^3 
-^44 — *4"222 ; 



7^(223) _ 2c)3 

J 44 — ''4"223 I 



^(233)_ 2^3 
J 44 — *4'-'233 ' 



(333) _ 2q3 



44 



^4^333 



(B.4) 



(B.5) 



One can show that the infinite-dimensional algebra generated by 43 generators R, B, T possesses 
infinitely many common invariant subspaces: it leaves invariant the space Vn for any n ^ Af and 
therefore preserves the flag 

Let us introduce an auxiliary generator 



jO = Sidi + 252^2 + 2S393 + 3S454 

Then one can define six 'raising' generators 

Ji+ = Sl J" , J3t-2 = 53^2 J° , 

•^22 -4 ~ S294 J , J23 _4 — 828-^843 , 



Jt _2 — 54^2 , 



J1 



81843"" . 



Finally, we get the infinite-dimensional algebra generated by 49 generators (B.2) 
called by definition f'^^\ 



(B.6) 
(B.6) and it is 



These raising generators (B.6) determine a highest-weight vector if n is a non-negative integer 
numb er. It l eads to a finite-di mens ional representation. It can be demonstrated that the operators 
(B.2) - (B^) leave the space (BJ.) invariant at fixed n a nd ac t on it irreducibly. Hence, as stated 
by the Burnside theorem fl^ , any ope rat or a cting on ( B.l ) allows a repres entat ion as a non- 
linear combination of the operators (B.2) - (B.6) plus an operator annihilating ( |B.l[ ) (annihilator). 
Therefore, the endomorphism of the space (B.l) is given by the infinite-dimensional algebra /^^^ 
gene rate d by 49 generators (B^) - (|B.(]| ). In turn, subalgebra of /^'^■' generat ed by the generators 
(B.2) - ( |B.5| ) possesses infinitely many finite-dimensional invariant subspac es ( B.l) at n = 0, 1, . . . 
and hence preserve the infinite non-classical flag of spaces of polynomials ( 3.17 ), T'q C Pi C 7^2 C 
• • • C P„ C . . . . 
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C Dual relations 



In this Appendix we present polynomial functions ga{^',P) which satisfy the eq.( 4.23| ) at ak = 
5fc(sin2(/3zfc/2)//32). These functions can be u sed as variables providing an algebraic form for the 
Hamiltonian in the z-representation ( 4.17 ). Compare these formulas with much simpler form 
of the functions — fa{cr; (3) from eg. ( 4. 11 ) which give an algebraic form to the Hamiltonian in 
x-representation. 

2 16 32 

93{<j;P) - -8<T3 + ^^i<J2 - + 4/5'(4a2<T? + ijt~ i2al + 120^1^3) 
o o io 



16 

g/3''(6(74<7i + 2ct3<72 + crlaz) + — /36(2ct4<T2 + <7i<74 + 2al) 



128 8- - , 128 10-2 
<^A<^3 + —P , 

2 4 12 

54(ct; P) = 16<T4 - 3 '5'i'^2 + + j^^t + (o-i(73 - 24d-i(74 - 4(730-2) 



^/3''( - (72^4 + I6.T4.T2 + a'i) - y/3^CT4<T3 + y 



56 (ct; /?) = yo-i - ^'^i'^2 - -^o"! + 16cti(T4 - 64(T4(T2 + ^^2^1 

+ i/3^(l92(T4CT3 + 192ct4(T2CTi - 48ct4(Ti - 160-30-2 + 8(T3(T20-i - (T3d-i) 
o 

2 

+ -/?■*(- 192ct^ - 9604CT3CT1 - 80d-4d-| + 160-40-20^ + o-4<Tf + 80-^0-2 
16 

- 2alaf) + —P^{72alai + 60040-302 - 604O30-1 - o|) 

+ |/3«( - 160^02 + ala! - 5a4a32) + ^P'Oaja, - ^/?^2o| . (C.7) 



D First eigenfunctions of the rational and trigonometric F4 
models 

In this Appendix we present explicit expressions for the first eigenfunctions of F4 models at n = 
0,1,2. 

I. Rational model. 
• n = 
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$0= 1, 
Eo = 0. 

• n= 1 

2 

*i = si (6At + 6t^ + l), 

(J 

£i = - 4a; . 

• n = 2 

$W=s? - -(4m + 4z/ + 1)si + -^(4;U + 4i/+1)(6m + 6i/+1) , 
E^^^ = - 8a; , 

^^2^ =S3 + ^(2/z + 4i/ + 1)s? - ^(2/x + 4z/ + l)(4/x + 4z/ + l)si 
+ ^(2/i + 4i/+l)(6M + 6z/+l)(4;U + 4z/+l) , 
^ = - 12a; , 

= S4 + -(3j/+1)s3 + ^(3:/+l)(2/x + 4:/ + l)s? 

- ^(3iv + l)(2;u + 4i/ + l)(4/i + 4i/ + l)si 

+ ^(3z^ + l)(2/x + 4i/ + 1)(6m + 6i/ + l)(4/i + 4z/ + 1) , 
) = - 16a; 



II. Trigonometric F4, model. 

• n = 

$0 = 1, 
Eo = 0. 

• n = 1 

, 6z/ + 6/i+l 

(5jv + 6^ + 1) 
iJi = 4(5z/ + 6m + l)/?^ . 

• n = 2 

(1)^ 3 (4i/ + 4M+l) 4 3 (4i. + 4M+l)(6i. + 6/i + l) 6 
2 8 (:^ + 4/z + l) 16 (z. + 4/i + l)(3z/ + 5/x + l) 

i;^^^ =8 (3 1/ + 5 /X + l)/?^ . 



32 




3 (4z/ + 4M+l)(4t/ + 2/x+l) 
16 (2t/ + 3^ + l) 



~ 16 (2j/ + 3At+ l)(3i^ + 4/x+ 1) 

12(3i/ + 4/u + 1)/?^ . 



_ 9 (3i/+l)(4z/ + 4M+l)(4t/ + 2M+l) 6 
64 (5i/ + 6/i + 3)(2z/ + ^ + l) ^ 

9 (3 z/ + 1) (4 + 4 + 1) (4 + 2 /X + 1) (6 + 6 /X + 1) 

128 (5 z/ + 6 M + 3) (2 1/ + /i + 1) (5 1/ + 6 M + 2) 

8(5i/ + 6m + 2)/J2 . 



3 (3t^+l) 
8(2!/ + /i + l) 



r'p3 + |:(z^ + l)/3-V 



33 



References 



[1] M. A. Olshanetsky and A. M. Perelomov, "Quantum completely integrable systems 
connected with semi-simple Lie algebras" , 
Lett. Math. Phys. 2 (1977) 7-13 

[2] M. A. Olshanetsky and A. M. Perelomov, "Quantum integrable systems related to Lie 
algebras" , 

Phys. Rep. 94 (1983) 313 

[3] F. Calogero, "Solution of a three-body problem in one dimension", J. Math. Phys. 10 (1969) 
2191-2196; 

F. Calogero, "Ground state of a one-dimensional A^-body problem", J. Math. Phys. 10 (1969) 
2197-2200; 

F. Calogero, "Solution of the one-dimensional A^-body problem with quadratic and/or 
inversely quadratic pair potentials", J. Math. Phys. 12 (1971) 419-436 

[4] A. V. Turbiner, "Lie algebras and linear operators with invariant subspace" , in Lie algebras, 
cohomologies and new findings in quantum mechanics (N. Kamran and P. J. Olver, eds.), 



AMS Contemporary Mathematics, vol. 160, pp. 263-310, 1994; f unct-an/9301001 



"Lie-algebras and Quasi-exactly-solvable Differential Equations", in CRC Handbook of Lie 
Group Analysis of Differential Equations, Vol.3: New Trends in Theoretical Developments 
and Computational Methods, Chapter 12, CRC Press (N. Ibragimov, ed.), pp. 331-366, 1995 
hep-th/9409068 

[5] W. Riihl and A. V. Turbiner, "Exact solvability of the Calogero and Sutherland models" , 
Mod. Phys. Lett. AlO (1995) 2213-2222 
liep-th/9506105 

[6] B. Sutherland, "Exact results for a quantum many-body problem in one dimension", Phys. 
Rev. A4 (1971) 2019-2021; 

B. Sutherland, "Exact results for a quantum many-body problem in one dimension, 11", 
Phys. Rev. A5 (1972) 1372-1376 

[7] L. Brink, A. Turbiner and N. Wyllard, "Hidden Algebras of the (super) Calogero and 
Sutherland models" , 
Journ. Math. Phy s.39 (1998) 1285-1315 
hep-th/9705219 

[8] M. Rosenbaum, A. Turbiner and A. Capella, "Solvability of the G2 integrable system" , 
Intern.Journ.Mod.P hys. A13, (1998) 3885-3904 



solv-int/9707005 



[9] A. Turbiner, "Hidden Algebra of Three-Body Integrable Systems" , 
Mod. Phys. Lett. A13 (1998) 1473-1483 



solv-int/9805003 



[10] O. Haschke and W. Ruehl, "Is it possible to construct exactly solvable models ?" 
Lect. Notes Phys. 539 (2000) 118-140 
hep-th/9809152 



34 



[11] N. Bourbaki, in Groups et Algebras de Lie (Hermann, Paris, 1968) Chaps. IV- VI, (V-5-4, 
prop. 5) 

[12] D. Bernard, V. Pasquier and D. Serban, "Exact solution of long-range interacting spin 
chains with boundaries," 
Europhys. Lett. 3 (1995) 301-305 
hep-tli/9501044 

[13] J. Wolfes, "On the three-body linear problem with three-body interaction," 
J.Math.Phys. 15 (1974) 1420-1424 

[14] A. Minzoni, M. Rosenbaum and A. Turbiner, "Quasi-Exactly-Solvable Many-Body 
Problems" , 

Mod. Phys. Lett All (1996) 1977-1984 
liep-th/9606092 

[15] A.V. Turbiner, "Quasi-Exactly-Solvable Problems and the SL(2,R) Group", 
Comm.Math.Phys. 118, 467-474 (1988) 

[16] X. Hou, M.A. Shifman, "A quasi-exactly-solvable iV-body problem with the sl{N + 1) 
algebraic structure" , 

Int.Journ.Mod.Ph ys. A14 (1999) 2993-3004 
liep-th/9812157 

[17] O. Haschke and W. Ruehl, "The construction of trigonometric invariants for Weyl groups 
and the derivation of corresponding exactly solvable Sutherland models ?" , 
Mod. Phys. Lett. A 14 (1999) 937-949 
inath-ph/9904002 

[18] S.P. Khastgir, A.J. Pocklington, R. Sasaki, "Quantum Calogero-Moser Models: Integrability 
for all Root Systems" , 
J. Phys. A33 (200 0) 9033-9064 
liep-th/0005277 

[19] S. Lang, "Algebra", Addison- Wesley Series in Mathematics, 

Addison-Wesley Publishing Company, Reading, Massachusets, 1965 



35 



